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ABSTRACT 

NGC 4473 is a so-called double sigma {2a) galaxy, i.e. a galaxy with rare, double 
peaks in its 2D stellar velocity dispersion. Here, we present the globular cluster (GC) 
kinematics in NGC 4473 out to ^ 10 Re (effective radii) using data from combined 
HST/ACS and Subaru/Suprime-Cam imaging and Keck/DEIMOS spectroscopy. 

We find that the 2cr nature of NGC 4473 persists up to 3 Re, though it becomes 
misaligned to the photometric major axis. We also observe a significant offset between 
the stellar and GC rotation amplitudes. This offset can be understood as a co-addition 
of counter-rotating stars producing little net stellar rotation. We identify a sharp radial 
transition in the GC kinematics at ^ 4 Rg suggesting a well defined kinematically 
distinct halo. In the inner region (< 4 Rg), the blue GCs rotate along the photometric 
major axis, but in an opposite direction to the galaxy stars and red GCs. In the outer 
region (> 4 Rg), the red GCs rotate in an opposite direction compared to the inner 
region red GCs, along the photometric major axis, while the blue GCs rotate along 
an axis intermediate between the major and minor photometric axes. We also find a 
kinematically distinct population of very red GCs in the inner region with elevated 
rotation amplitude and velocity dispersion. 

The multiple kinematic components in NGC 4473 highlight the complex formation 
and evolutionary history of this 2cr galaxy, as well as a distinct transition between the 
inner and outer components. 

Key words: galaxies: star clusters - galaxies: evolution - galaxies: kinematics and 
dynamics - globular clusters 


1 INTRODUCTION 

Kinematically distinct components (KDCs) in early-type 
galaxies (ETGs) provide clues about their formation and 
evolution. The KDCs can be viewed as fossils of the accre¬ 
tion and merger processes that built up present day galax¬ 
ies in the hierarchical formation framework. Early studies of 
the kinematics of ETGs were limited to the central regions, 
where rotation, sometimes along a misaligned or twisted ro¬ 
tation axis, as well as kinematically distinct components 
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were observed ( 

Davies et al. 2001 Emsellem et al. 2004 

Krajnovic et al. 

2011). More recently, the probed radii have 


been extended into the halo region, where some of these fea¬ 
tures have been shown to persist, while in others, changes 


are seen at large radii ( 

Arnold et al. 2014 Raskutti et al. 

2014 Foster et al.||20I5 

). Stellar light is, however, faint in 


the halo region, hence the need for bright kinematic tracers 
like planetary nebulae and globular clusters (GCs). 

Galaxies, apart from forming GCs “in-situ”, are ex¬ 
pected to have acquired some GCs formed “ex-situ” through 
galaxy mergers in the two-phase galaxy formation model 
(e.g. Oser et al.|2QI0 ). Almost all ETGs, studied with deep 
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enough photometry, have been shown to have at least a bi- 


modal GC colour distribution (e.g. Zepf & Ashman 1993 
Ostrov et al.|1993 ). This colour bimodality usually points at 


an underlying bimodality in metallicity ( Usher et al.|2012 ), 
which has been linked with the hierarchical merging his¬ 


tory of the host galaxies (Tonini 2013). It is possible, by 


observing correlations in position-velocity parameter space 


of GCs (Romanowsky et ah 2012 Blom et ah 2014 Fos¬ 


ter et ah 2014), to unearth relics of the assembly history 


of the host galaxy. This is because in the galaxy outskirts, 
dynamical timescales are longer - on the order of ~ 1 Gyr 
(Goccato et al. 2013), hence a “memory” of the infalling 


galaxies’ orbital properties is expected to be retained (John- 
ston et aL]|1996 Helmi White|p~999 ). GGs therefore con¬ 


tain fossil records of the chemical and dynamical processes 
that shaped the structure of present-day galaxies. Similari¬ 
ties between the kinematics (amplitude of rotation velocity, 
velocity dispersion, direction of rotation, etc.) of stars and 
GGs in galaxies ( Foster et al.|20lT||Pota et al.|2013 ) can be 
used to further constrain their likely assembly history. 

Double sigma {2a) galaxies ( Krajnovic et al.||2011 ) are 
examples of multi-spin galaxies ( Rubin||1994 ), easily identi¬ 
fied by observing their stellar 2D velocity dispersion maps. 
These maps uniquely show a pair of symmetric, off-center, 
significant bumps {2a peaks) aligned along the photometric 
major axis ( Bois et al.pOllj |Krajnovic et ah' 2011 


Foster 


et al.|20i3 ). This feature has been linked to the presence of 


a pair of extended, counter-rotating disc-like stellar struc¬ 
tures at the center of the galaxy. 2a galaxies are rare - ~4% 
in the ATLAS^^ survey ( Krajnovic et al.|[2011 ) - probably 
because the conditions required to produce them are uncom¬ 
mon, as shown in idealised binary merger simulations (e.g 
Jesseit et al. 1|2007| |Bois et al.||mll l. They were formed in 


these simulations only when the merging spiral progenitors 
have similar masses (mass ratio 1:1 to 3:1) and are co-planar 
( Grocker et al.|2Q09 ). It should be noted that an alternative 
channel for forming galaxies with counter-rotating stellar 
components is accretion of cold gas along cosmological fila- 


ment s ( Rix et al.|1992[[Thakar Ryden|19% Algorry et al. 


2014). Generally, the 2a galaxies from the ATLASsurvey 


are flattened (ellipticity, e > 0.4), viewed at high inclination 
{i > 70°) and have low to intermediate luminosities. 


NGG 4473 is the most massive known 2a galaxy (Kra¬ 


jnovic et al. 2011), though it can be described as an L* 


galaxy. It has an absolute magnitude of Mr = —23.8 and 
is photometrically undisturbed. It is a flattened (e = 0.43) 
elliptical galaxy at a distance of 15.2 Mpc in the Virgo Glus- 
ter, with effective radius Re of 27" and photometric position 
angle of 92.2° ( Brodie et al.|[20T4 and references therein). 
It has a galaxy recession velocity (Usys) of 2260 km . 
A Sersic function fit to its surface brightness shows that 
NGG 4473 is a “cuspy” galaxy with excess central light (Ko- 
[rmendy et al.|200'9 Dullo Graham|2013 ), suggesting a dis¬ 
sipative merger origin. It has fast rotation in its inner region 
(< 1 Re) and harbours a pair of counter-rotating stellar discs 
( Gappellari et al.[|2007 Krajnovic et al.||2011 ). Foster et al. 
(2013) studied the stellar kinematics out to ~ 3 Re and ob¬ 
served counter-rotation which extends beyond ~ 60^^ along 
the photometric major axis. They used the SKiMS (Stellar 
Kinematics with Multiple Slits) method ( Norris et al.|20(3^ 


Proctor et al. 2009) which extends the stellar kinematics 


lar kinematic components in the galaxy outskirts, rotating 
along both the photometric major and minor axes. They 
concluded therefore that NGG 4473 is triaxial, with a kine¬ 
matically distinct halo (KDH). 

Here, we study the kinematics of the GG system of NGG 
4473 usi^ photometric and spectroscopic data from the 
SLUGG^(SAGES Legacy Unifying Globulars and Galax- 
ieS) survey ( Brodie et al.|[20T4 ) out to large galactocentric 
radii (i.e. ~ 10 Re). While Foster et al. (2013) were able 
to identify a KDH, they could only probe the inner edge of 
this region. GGs are better suited to probe this KDH prop¬ 
erly. We therefore study the radial profile of the GG system 
kinematics, with particular interest in any sharp kinematic 
transition(s). These transitions can be used to understand 
the nature of the last major merger the galaxy experienced, 
as proposed by Hoffman et al. (2010). We also investigate 


the radial extent of the 2a velocity dispersion feature, using 
GGs. We study the kinematics of the blue and red GG sub¬ 
populations separately: differences in their kinematics, be¬ 
sides providing additional evidence for GG colour bimodal¬ 
ity, also provide clues about the progenitor galaxies (|Bekki 
et al.||2005| ). 

This paper is structured as follows: In Section we 
present the photometric and spectroscopic data used in this 
study. In Sectionj^ we focus on the photmetric analysis, pro¬ 
duce 2D mean velocity and velocity dispersion maps, as well 
as statistically determine the significance of the 2a feature. 
We also produce ID kinematic radial and colour profiles. 
Section ends with a discussion of the line of sight veloc¬ 
ity distributions of the GG system. In Section we briefly 
summarise models from the literature which form 2a galax¬ 
ies. In Section we discuss our results, relating them to 
model predictions from the literature. Finally, in Section 
we briefly summarise the paper. 


2 OBSERVATIONS AND DATA REDUCTION 

NGG 4473 was observed with Suprime-Gam on the Subaru 
telescope on the night of 2010 November 4. The total ex¬ 
posure times were 688, 270 and 450 seconds with average 
seeing of 0765, 0765, 0767 in the g, r, and i bands, respec¬ 
tively. The raw images were reduced using the SDFRED2 re¬ 
duction pipeline ( Ouchi et al.|2004 ). To detect GGs, a model 
of the galaxy light using the IRAF/ellipse task was first 
subtracted from the reduced images and the IRAF/DAOFIND 
task was then used to find bright and compact objects. At 
a distance of 15.2 Mpc, GGs in NGG 4473 are unresolved 
and appear as point sources. The objects detected from dif¬ 
ferent bands were matched and aperture photometry was 
carried out with the IRAF/phot task to remove extended 
sources using the method described in Spitler et al. (2008). 


Finally, we corrected our photometry for reddening using 


the dust extinction maps from Schlegel et al. (1998). We 


used a {g — r) vs {r — i) colour-colour plot to identify GG 
candidates. We supplement our photometric catalogue with 
g and z band photometry from the Hubble Space Telescope 
{HST) Advanced Gamera for Surveys (AGS) for NGG 4473 
(see Strader et al. [20(36 and Spitler et al.|2006 for details of 


to relatively large radii. They also identified multiple stel- 
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the photometric reduction). We used the colour transforma¬ 
tion equation from [Usher et ah (2012) to convert the HST 
z band photometry into the i band. Our final photometric 
catalogue consists of 1097 GC candidates. 

We used the DEep Imaging Multi-Object Spectrograph 
(DEIMOS) on the Keck II telescope to obtain spectra for ob¬ 
jects that are probable GCs from our photometric catalogue. 
The spectroscopic observations were taken on the nights of 
2011 March 30, 2012 Eebruary 19 and 2012 Eebruary 20. See¬ 
ing was between 0.7" - 0.95" and four DEIMOS masks were 
observed during the campaign. We used the usual SLUGGS 
setup described in Pota et al. (2013) and integrated for 2 
hours per mask, though one of the masks was observed for 
2.75 hours. We reduced our raw spectra using the DEEP 
spec2d pipeline ( Gooper et al.j2012 ) in IDL. We determined 
the heliocentric radial velocities by cross-correlating each 
spectrum with 13 stellar templates, obtained with the same 
DEIMOS instrumental setup, using the IRAF/RV. FXCOR task. 
Eor each object, the radial velocity is the average of the mea¬ 
sured radial velocities from FXCOR. The uncertainty of each 
radial velocity is estimated by adding in quadrature the er¬ 
ror output from FXCOR to the standard deviation among the 
templates, which is an estimate of the systematics. 

The peaks of the Galcium Triplet (GaT) lines (at 8498, 
8542 and 8662 A) and the Ra (6563 A) line were used to 
identify GGs associated with NGG 4473 based on their radial 
velocities. Eor secure GG confirmation, we require at least 
two clearly visible GaT lines (usually 8542 and 8662 A) and 
the Ra when probed, as well as GG-status consensus by at 
least two members of the team. Eour spectroscopic GGs had 
repeated measurements from different masks, all in good 
agreement and consistent within the uncertainties associ¬ 
ated with our observations. The final measurements used 
for these objects (and recorded in our final spectroscopic 
object catalogue) are the weighted average of the individual 
measurements. These repeated observations show that the 
errors in our radial velocity measurements are <15 km s“^. 
Our final spectroscopic catalogue (see Table contains 105 
unique GGs with redshift-corrected radial velocities, along¬ 
side Galactic stars and background galaxies. 

In Eig. we show the radial velocity distribution of 
GGs as a function of galactocentric radius. We have used 
the friendless algorithm of Merrett et al. (2003) to identify 
possible outliers from our spectroscopic GG catalogue. We 
implemented the algorithm to identify objects with radial 
velocities outside the 3a envelope of their 20 nearest neigh¬ 
bours. The algorithm returned no outliers. We note here the 
apparent asymmetry in the velocity distribution in the in¬ 
ner 2', especially in the red GGs. This can be viewed as a 
signature of rotation, an issue we explore in details later in 
this paper. 



Figure 1. Globular cluster (GC) radial velocity distribution as 
a function of galactocentric radius. The solid line and dashed 
curves show the recession velocity of NGC 4473 and the 3a enve¬ 
lope of GC velocities, respectively. Blue circles and red stars are 
spectroscopically confirmed blue and red GCs, respectively. 
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(s - t) 

Figure 2. Globular cluster (GC) {g — i) colour distribution. 
Top panel: colour magnitude diagram (CMD) of GC candidates 
from HST and Subaru photometric catalogues. Black, open and 
filled circles are HST and Subaru GC candidates, respectively. 
The vertical line (in both panels) shows the fiducial colour cut 
{9~i) = 0-93 mag used to separate blue and red GCs in our spec¬ 
troscopically confirmed sample. Mean photometric uncertainties 
are shown on the right as error bars. Bottom panel: Black his¬ 
togram is the distribution of GC colour from combined HST and 
Subaru photometric catalogues, while the red histogram is the 
colour distribution of spectroscopically confirmed GCs. Gaussian 
fits from GMM are over-plotted on the histogram. The GC colour 
distribution is bimodal at a 99.99% significance level. 


3 ANALYSIS 

3.1 Photometric Analysis 

The final colour-magnitude diagram (CMD) of all unique 
GC candidates from the combined HST and Subaru data, 
brighter than i < 23.5 mag - one magnitude fainter than the 
GC turnover magnitude, Mtomg = —8.4 mag from jVille-j 
- is shown in the upper panel of Fig. 


gas et al. (2010 


We use an upper magnitude cut of i > 19.9 mag (using the 
integrated magnitude of cj-Cen, Mi ^ —11.0 mag from [Pot aj 
et al.|2013 and distance modulus from Villegas et al.|2010| ). 


to isolate likely ultracompact dwarf (UCD) candidates from 
our sample, since they may have different chemo-dynamical 
properties (see Strader et al.|20rT ). While this cut excludes 
the two brightest GCs in our sample (with z-band magni- 
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ID 

RA 

Dec 

U 

AU 

9 

Ag 

r 

Ar 

i 

Al 


[Degree] 

[Degree] 

[km s“^] 

[km s“^] 

[mag] 

[mag] 

[mag] 

[mag] 

[mag] 

[mag] 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

NGC4473_GC1 

187.45086 

13.42309 

2112 

14 

22.492 

0.016 

21.905 

0.016 

21.640 

0.017 

NGC4473_GC2 

187.45636 

13.42486 

2390 

12 

22.024 

0.013 

21.387 

0.013 

21.102 

0.014 

NGC4473_GC3 

187.46807 

13.42744 

2524 

20 

21.715 

0.009 

20.953 

0.009 

20.569 

0.008 

NGC4473_GC4 

187.46594 

13.42724 

2327 

12 

23.364 

0.022 

22.679 

0.022 

22.309 

0.024 

NGC4473_GC5 

187.46444 

13.43094 

2279 

12 

22.251 

0.011 

22.682 

0.012 

21.442 

0.013 

NGC4473^tarl 

187.47899 

13.52667 

-206 

12 

22.992 

0.011 

22.391 

0.011 

22.196 

0.011 

NGC4473^tar2 

187.38470 

13.53626 

-148 

23 

22.609 

0.009 

22.095 

0.010 

21.915 

0.010 

NGC4473_gall 

187.42701 

13.48026 

z = 0.3 

— 

23.067 

0.011 

22.464 

0.011 

22.284 

0.012 

NGC4473_gal2 

187.52558 

13.36814 



22.595 

0.010 

22.107 

0.010 

21.863 

0.010 


Table 1. Spectroscopically confirmed globular clusters, stars and galaxies. Notes Column (1): Object identifier, written as galaxy name 
and object type (globular clusters - GC, stars and galaxies). Columns (2) and (3): Object position (RA and DEC, respectively) in 
degrees (J2000.0). Columns (4) and (5): Measured heliocentric radial velocities and uncertainties, respectively. Columns (6)-(ll): Subaru 
gri photometry and corresponding uncertainties. The photometry has been corrected for Galactic extinction. The full version of this 
table is available in the online version. 


tudes of 19.8 and 19.85), we however chose to retain them 
after cross-matching them with the published catalogue of 
GCs in NGC 4473 from [Jordan et ah ( |2009| , where they 
have sizes consistent with being GGs. 

Using the Gaussian Mixture Modelling (GMM) code by 
Muratov & Gnedin (2010), we determine that the GG colour 


distribution of NGG 4473 is best described by a bimodal, 
heteroscedastic distribution at a 99.99% significance level. 
This is shown in the histogram in the bottom panel of Fig.|^ 
where the modes are clearly well separated. The (g — i) peaks 
for the Gaussian modes are 0.78 and 1.06 with dispersions 
of 0.08 and 0.09 respectively. We group our GG candidates 
into blue and red GG subpopulations, using a fiducial colour 
split of (g — i) = 0.93 mag, and find that the fractions of 
blue and red GGs in our photometric sample are 0.58 and 
0.42 respectively. 


3.2 Mean velocity and velocity dispersion 2D 
maps 


We use the method of Goccato et al. (2009) and Pota et al. 


(2013) to construct mean radial velocity and velocity dis¬ 


persion 2D maps for the spectroscopic GG sample. We con¬ 
struct an equally spaced nx?n=100xl00 grid in position 
space and at every grid point, we use the iV— nearest GGs 
to compute the interpolated mean velocity, v(i,j) and ve¬ 
locity dispersion, cr(z,j) (see Eqns. and |^, weighting the 
measured radial velocities Vk{a,S) by the inverse square of 
their separations from the grid point, Wk- The optimal num¬ 
ber of the nearest neighbours, A" = 10, used was determined 
as the square-root of the GG sample size, as suggested by 


Pinkney et al. (1996). We however varied this by 20% to 


test the robustness of our maps (i.e. A = 8,12) and found 
no significant difference in the map features. It should be 


noted that the maps we show here are mostly illustrative, 
hence we only quote qualitative trends from them. 




Y,Vk{a,5)/wk^ 

N 


( 1 ) 




J2Vk(a,Sf/wk^ 

k=l 




1/2 


( 2 ) 


and 


In Eqns. and (i^j) and (a, (5) are grid 
sky coordinates respectively, Av(i,j) is obtained as in 
E qn. but using th e estimated uncertainties and w = 
y /(a — i)^ (S — j)2. The 2D maps shown in Eigs. and 
1^ have additionally been smoothed with Gaussian kernels 
that vary linearly with radius to recover the global kinemat¬ 
ics. We produce mean velocity and velocity dispersion 

2D maps for combined SAURON ( Emsellem et al.|2004 ) 
and SKiMS data ( Poster et al.||2Q13 ) for comparison. Even 
though the datasets have different spatial extents and sam¬ 
pling, the maps agree well. In the left panel of Eig. 
where we have combined SAURON and SKiMS data, NGG 
4473 shows an axisymmetric, disc-like structure that ex¬ 
tends to the edges of the SAURON map. The dumb-bell 
feature in the velocity dispersion map (right panel of Eig. 
1^ is typical of 2a galaxies and the peaks are separated by 
~3 Re- The kinematic position angle (PAkin), defined as 
the angle between the direction of maximum rotation and 
North, is aligned with the photometric major axis (92.2°). 
There is significant minor axis rotation beyond 50^^ as well 
as counter-rotation along the photometric major axis (as 
reported by Poster et al. 2013). Note that rotation along 


an axis e.g. major axis, implies that the velocity signature 
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Figure 3. 2D smoothed mean velocity (left panel) and velocity dispersion (right panel) maps for galaxy stars in NGC 4473, using data 
from SAURON ( [Emsellem et al.]|2004| ) and SKiMS ( [Foster et al.|20T^ . In the left panels rapid disc-like rotation and counter-rotation 
beyond the footprint of the SAURON data (i.e. < 1 Re) can be seen along the photometric major axis. In the right panel, the 2(7 feature 
is prominent along the major axis. Stellar isophotes at I and 2 Re have been over-plotted in both panels. 




Figure 4. 2D smoothed mean velocity (left panels) and velocity dispersion (right panels) maps for globular clusters (GCs) in NGC 4473. 
In the top, middle and bottom panels, we show all the GCs, the blue GCs and red GCs, respectively, using the method described in Sec. 
|3.2| The ellipses are 2 and 7 Re stellar isophotes and the circles correspond to observed GC positions. In the top right panel, there are 
signs of the 2a feature in the velocity dispersion map for all GCs out to ~ 3 Re, although misaligned to the photometric major axis. In 
the middle left panel, the blue GCs rotate intermediate between the photometric major and minor axes in the galaxy outskirts, while in 
the bottom panel, the red GCs counter-rotate along the photometric major axis, similar to the galaxy stars as shown in Fig. 
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would be seen around the perpendicular axis i.e. the minor 
axis. In Fig. ^ the red GCs dominate the rotation of the 
GC system in the inner ~7 Re region, resulting in 2D maps 
akin to that of the galaxy stars, with counter-rotation ev¬ 
ident along the photometric major axis (bottom left panel 
of Fig. [^. The mean velocity map of the blue GGs (middle 
left panel) shows a non-regular rotation pattern, especially 
in the outer region, with rotation along an axis intermediate 
to the photometric major and minor axes. The velocity dis¬ 
persion maps for the blue and red GGs are similar, with the 
red GGs having higher velocity dispersion within the inner 
~7 Re region. However, we see a 2cr-like feature in the veloc¬ 
ity dispersion map when all the GGs are combined together 
(top right panel). These elevated, offset bumps in the 2D 
velocity dispersion map are however not aligned along the 
photometric major axis. 


3.3 Significance of the extended 2a feature 


To determine how far out the 2a feature extends in the 
galaxy stars and GGs, we assume that our measured ve¬ 
locities (and their associated errors), at all positions on the 
sky, are drawn from the same distribution (see Walker et al. 
2QQ6| . We then create mock catalogues (Wzm = 1000), using 
the same sky positions as in our real catalogue, but with ve¬ 
locities (and errors) drawn randomly, without replacement, 
from the measured values. This non-parametric approach 
ensures that each mock catalogue has the same velocity dis¬ 
tribution as the measured data, but destroys any correla¬ 
tion between position and velocities. We use the method 
in Sec. |3.2| to then make velocity dispersion maps for each 
mock catalogue, and note the maximum velocity dispersion, 
c^mock per mock catalogue. The significance of the 2a fea¬ 
ture at every grid point p(z, j) is then quantified as the ratio 
N{areai{i,j) > cTmock) /Ngim■ From Fig.the 2a feature is 
significant in both galaxy stars and GGs with p > 0.95. The 
symmetrical form of the 2a feature is observed in the GGs as 
well, within 3 Re, though it appears inclined at ~ 30° to 
the photometric major axis. This misalignment needs to be 
studied in more detail to understand its true nature. How¬ 
ever, from the mean velocity maps in Fig.|^ a large velocity 
difference is evident between the blue and red GGs in the 
spatial region corresponding to the 2a feature. Also, these 
velocity offset regions are misaligned, with the blue GGs ro¬ 
tating along an axis intermediate to the photometric major 
and minor axis, coincident with the axis of the misaligned 
2a feature. We therefore conclude that the blue GGs drive 
the observed misalignment in the 2a feature. 


3.4 Kinematic ID radial profiles 

To understand the one dimensional (ID) kinematic radial 
profiles of the GG system, we fit an inclined disc model to 
our discrete GG data ( Foster et al.||20lT| . We evaluate the 
circularized galactocentric radius, R using: 


R = 


^photA^^ T 


Qphot 


(3) 


where ^phot is the ratio of the minor to major photometric 
axis (^phot = 1 —e), and X and Y are the projected cartesian 



100 50 0 -50 -100 

ARA [arcsec] 



400 200 0 -200 -400 

ARA [arcsec] 


Figure 5. Contours of the 2a velocity dispersion feature in 
galaxy stars {top panel) and GCs {bottom panel). Contour levels 
are taken from the fraction of mock catalogues with a maximum 
velocity dispersion less than the estimate from the observed ve¬ 
locity dispersion. Blue and red CCs are overplotted in the bottom 
panel as blue circles and red stars, respectively. The ellipse in 
the top panel is the 2 Re stellar isophote. In the bottom panel, we 
show the 2 and 7 Re stellar isophotes. The 2a feature is significant 
in the CCs, along an axis inclined at 30° to the photometric 
major axis, out to ~ 3 Re- 


coordinates of individual GCs on the sky, relative to the 
galaxy photometric axes. The model we fit is : 


kmod.'i — V^ys i 


Hrc 


1 + 


f tan(PA^— 




(4) 


Foster et al. (2011), Blom et al. (2012) and Pota et al 


(2013), where we minimise the function: 


x"oc^ 


(Vi - Vmod,if 

(a2 + (AVi)^) 


+ ln(cr^ + (At7)^) 


(5) 


to obtain the best fit amplitude of rotation Kot, ve¬ 
locity dispersion a and kinematic position angle PAkin- In 
Eqns. 1^ and Vsys is the recession velocity of the galaxy, 
while Vi, AVi and PAi are the radial velocities, uncertainties 
in measured radial velocities and position angle for the 
GC, respectively. We do not fit for ^kin, the flattening due to 
rotation, since it is not easily constrained with sparse data. 
We therefore consider two cases: (i) where we fix gkin = ^phot 
and (ii) ^kin = 1- Both cases result in similar GC kinematics, 
hence we adopt profiles where gkin = ^phot, as done by |Foster| 
et al. (2013). We fit kinematic parameters by minimising the 
function in Eqn. in rolling bins of fixed binsize (A=20), 
starting with the innermost GCs. This choice of binsize en- 
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R/Re 



Figure 6. Kinematic radial profiles for blue and red globular clusters (GCs) and stars in NGC 4473. Top, middle and bottom panels 
show the rotation velocity (Krot), velocity dispersion (cr) and kinematic position angle (PAkin) from the kinematic fit, respectively. We 
show 68 percent confidence interval contours in all panels. Blue (solid), red (dashed), green and black contours are for blue GCs, red 
GCs, SKiMS and SAURON data, respectively. Note the offset in rotation amplitude between galaxy stars and GCs in the top panel. 
In the middle panel, the blue GCs have a flat velocity dispersion profile. In the bottom panel, the dashed and dot-dashed lines are the 
photometric major and minor axes, respectively. The galaxy stars have major and minor axis rotation while the red GCs at < 4 Re 
rotate in the same sense as the galaxy stars. There is significant counter-rotation along the major axis between the inner blue and red 
GCs up to ~ 4 Re. Beyond this, the red GCs begin to wrap around the North direction, rotating along the major axis, but in opposite 
direction to the galaxy stars. The blue GCs can be seen to rotate predominantly along the minor axis in the outer regions. 


sures that K-ot is not artificially inflatecQ especially when 
rotation is small. To obtain the 68 percent confidence in¬ 
terval on our kinematic parameters, we randomly sample 
1000 times (with replacement) from each bin and minimize 
the function in Eqn. The bin is continuously updated by 
deleting the innermost object and adding the next further 
out GC until we get to the outer radial boundary. 

In Fig. we show radial profiles of the kinematic pa¬ 
rameters for SAURON, SKiMS, blue and red GCs data for 


The amplitude of velocity rotation is artificially inflated when 


Vrot/o- < 0.55(Strader et al. 


20111. 


NGC 4473. The kinematic parameters for SAURON and 
SKiMS data were obtained using the method described in 


Foster et al. (2011). We fit for U-ot and PAkin together, and 


<j separately, keeping gkin = ^phot for consistency. We have 
applied an offset of 20 km s“^ to the SAURON velocity dis¬ 


persion as discussed in Arnold et al. (2014) and Foster et al. 


(2015). We constrained U-ot and a to vary between 0 and 
300 km s“^ while PAkin was allowed to vary freely in order 
to probe all possible kinematic components. 

The model we fit may not be sensitive to the radially 
superimposed components identified in the stellar kinemat¬ 


ics study of Foster et al. (2013), since it fits a single com¬ 


ponent function at every radius. To unravel the likely de- 
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Figure 7. Same as in Fig.|^ but for stars and all GCs in the inner 
< 4 i^e, with the purple lines now showing profiles for all the GCs, 
binned together in this region, regardless of their colour. We show 
the Ifj uncertainties for SAURON, SKiMS and all GCs. In the 
bottom panel, we do not show the uncertainties for clarity. The 
rotation amplitude, when the GCs are not separated by colour, 
is now significantly reduced with respect to that of the blue and 
red GCs, and now comparable to that of the stars. 


generacy in the kinematics, we therefore study kinematics 
in the blue and red GCs, separately. While contamination 
from the overlapping tails of the blue and red GC distribu¬ 
tions could introduce some uncertainty in our fit, we expect 
our fit to be robust enough to show dominant kinematic 
trends. As a further test of the robustness of our fit, we ob¬ 
tain kinematic profiles for the GC subpopulations, adopting 
only GCs with classification probability greater than 0.68. 
We recover profiles similar to those obtained when we ap¬ 
plied a straight colour cut, confirming that our fitting is 
indeed robust against the overlap. We also tested the ro¬ 
bustness of our fit against the choice of N^ the number of 
GCs per bin. We varied N by 20% (i.e. N — 16, 24) and 
found that our kinematic profiles are similar in all cases. 

As seen in Fig. the blue and red GCs show similar 
rotation amplitude profiles at all radii. However, there is 
a large and significant offset between their rotation ampli¬ 
tudes compared to that of the galaxy stars. This is strange 


and contrary to results from previous studies (e.g. Pota et al. 


2013) of other GC systems where the stellar rotation profile 


usually matches with that of the red GCs, though a sim¬ 


ilar discrepancy was observed in NGC 4494 (Foster et al. 


2011). We tentatively associate this offset with the counter¬ 


rotating stellar disks believed to be responsible for the 2cr 
feature - an issue we address later in Sec. [H Here we inves¬ 
tigate further by modelling the kinematics for all the GCs 
(i.e. combined blue and red GCs) in the region where we 
have very good azimuthal sampling and Kot is significantly 
high (< 4 i^e), and show the result in Fig. Hrot for all 
GCs is significantly reduced, and comparable to that of the 
galaxy stars, giving credence to our earlier hypothesis. 

The velocity dispersion profiles for both GC subpopu¬ 
lations are similar in the inner ^ A Re region. We see signs 
of decline in the velocity dispersion of the red GCs while 


the blue GCs maintain a flat profile. In the region of over¬ 
lap (see Fig. 0, the velocity dispersion profiles of the red 
GCs and galaxy stars are similar. Both GC subpopulations, 
however, differ in their kinematic position angle profiles. In 
the inner 4 i?e, the blue and red GCs counter-rotate along 
the photometric major axis, with the red GCs rotating in 
the same sense as the inner galaxy stars. The red GCs in 
the outer region (> 4 Re) also rotate along the photometric 
major axis but counter-rotate with respect to the inner red 
GCs. In the outer region, the blue GCs show rotation along 
the minor axis, with multiple components. Hence, we a see 
counter-rotation along the photometric major axis as well 
as a minor axis rotation in the GC kinematics, similar to 


the result of Foster et al. (2013) for the stellar component. 


In the study of GC kinematics, it is possible to over¬ 
estimate the rotation amplitude ( Sharpies et al.||1998 Ro- 


manowsky et al. 2009 Strader et al. 2011| Zhang et al. 2015), 
and we therefore perform statistical tests to ascertain that 
our rotation profile is real. We quantify possible bias using 
GCs in the inner (< 4 Re) radial bin by generating artificial 
discrete velocities at the GC position angles using our kine¬ 
matic model. In the model, we fix a and PAkin to the best-fit 
values for this region and vary Kot over the range of possible 
values (from 0-300 km s“^). Vrot retrieved from the artificial 
dataset reveals that our kinematic fitting only suffers from 
bias at extreme rotation amplitudes, i.e. we tend to over¬ 
estimate Hrot at low intrinsic rotation amplitudes (e.g. at 
30 km s“^, we over-estimate Id-ot by 30 km s“^) and under¬ 
estimate at very high rotation amplitudes (e.g. at 270 km 


values for the blue GCs, the bias in Hrot is negligible while, 
for the red GCs, it is ~10 km s“^. 

We also created artificial datasets by randomly shuf¬ 
fling the position angles of the GCs in the inner radial bin, 
and then fit for kinematic parameters. We constructed 1000 
such datasets. The random shuffling serves to destroy any 
intrinsic correlation between GC position angle that could 
be driving the rotation amplitude. From the fit, we deter¬ 
mined the probabilities that the GC rotation is non-zero for 
the blue and red GC in the inner radial bin are 95 and 89 
%, respectively. 

We show the rotation dominance parameter {Vrot/cr) 
in the top panel of Fig. This parameter quantifies how 
“discy” a system is, such that when Hrot/cr > 1, the system 
is described as rotation-dominated. Though the blue and red 
GCs show high rotation within the inner ~ 4 Pe, the GC 
system is overall not rotation-dominated. The Vrot/cr profile 
has a negative gradient, at least within the inner ~ 4 Pe and 
we note the spike in the profile for the blue GCs at ~ 6 Pe- 
This may be a substructure signature. 

We also show the Hrms profile for blue and red GCs in 
Fig.|^ obtained using: 

1 ^ 

= N - ^sys)" - (AVi)" (6) 


^ rms 


and evaluating the uncertainty on our estimated Vrms using 
the relation from Danese et al. (1980). The Hrms is equiv¬ 


alent to the velocity dispersion in the absence of projected 
rotation and it is a measure of the total specific kinetic en¬ 
ergy of the galaxy. Overall, the profile for red GCs declines 
with radius, unlike the blue GCs which have a flat radial 
profile, with two bumps at ~ 3 and ~ 6 Pe- Since the blue 
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Figure 8. Rotation dominance, Rrot/cr {top panel) and Rrms 
{bottom panel) profiles for blue GCs and red GCs. Golour scheme 
is the same as in Fig. Below the dashed horizontal line, Vrot 
would be artificially inflated for our binsize choice (see Sec. |3.4| 
for details). The V-rot/cr profiles suggest that both blue and red 
GGs have high rotation. The TTms profile for the red GGs declines 
more rapidly than that of the blue GGs, which stays relatively flat 
at all radii. 


and red GCs independently trace the same galaxy potential, 
differences in their Vrms profiles are linked to their different 
orbital anisotropies and spatial distributions. 
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Figure 9. Globular cluster kinematics as a function of {g — i) 
colour in the inner < 4 Re (black-dashed contours) and outer 
> 4 Re (green-solid contours) regions. Dashed and dot-dashed 
lines are the same as Fig. and contours show the 68 percent 
confidence interval. The black vertical line is the colour cut, {g — 
i) = 0.93 mag used in our analysis. GCs in the inner bin show 
counter-rotation along the photometric major axis while the outer 
bin is dominated by minor axis rotation. The inner region has a 
population of very red GCs i.e. {g — i) >1.05 mag with high 
velocity dispersion that are rotating along the major axis. 


3.5 Kinematics as a function of GC colour 

The sharp twist in P^kin for both blue and red GCs at 
~ 4 Re (Sec. |3.4| ) suggests that this could be a radial bound¬ 
ary between the inner and outer GC subpopulations. We 
therefore use this fiducial radial limit to define inner and 
outer GC subsamples and study their kinematics as a func¬ 
tion of colour. For both cases, we use the model and method 
in Sec. |3.4[ binning by GC colour and fitting for the kine¬ 
matic parameters while fixing gkin = ^phot- We show the 
results in Fig. and note the similar Vrot profiles for GCs 
with {g — i) < 1.05 mag. For GC colours redder than this 
in the inner radial bin (i.e. < 4 Re), K-ot and a are signifi¬ 
cantly higher than the average for the entire GC system and 
P^kin is aligned along the photometric major axis. Again, 
the inner blue and red GCs counter-rotate along the photo¬ 
metric major axis. In the outer region, for the blue GCs, the 
dominant rotation mode is along the minor axis. 


3.6 Line—of—sight velocity distribution 


To better understand the variations of GC kinematics with 
radius and colour, as well as the sharp transition at ~ 4 Re, 
we construct line-of-sight velocity distributions (LOSVDs) 
for our spectroscopic GC dataset. We consider three sam¬ 
ples: all GCs in a single bin, GCs within ~ 4 Re in a single 
bin and GCs beyond ~ 4 Re from the galaxy centre in a 
single bin. We make LOSVDs for the blue and red GC sub¬ 
populations using the bins defined above and present the 
results in Fig. The LOSVDs have been smoothed with 
optimal bandwidth Gaussian kernels ( Silverman||1986 ) and 
normalised for easy comparison. 

From the top panel, all our three samples have LOSVDs 
that are approximately Gaussian with peaks consistent with 
the galaxy recession velocity. This can be interpreted as 


overall dynamic equilibrium. The middle panel, however, 
shows that the red GCs in the inner radial bin have an 
asymmetric LOSVD with a significantly redshifted veloc¬ 
ity peak and an excess of low-velocity GCs. The LOSVD 
is such that the prograde wing has a steeper gradient than 
the retrograde wing. In the outer bin, the LOSVD has a 
flattened shape, which might be due to tangentially-biased 
orbits. This suggests that the inner red GCs are dynamically 
different from those in the outer radial bin. In the bottom 
panel, the “peaked” shape of the LOSVDs suggest that the 
blue GCs, especially in the outer bin, are likely to be on 
radially-biased orbits. 


4 MODEL PREDICTIONS FOR 2a GALAXIES 


Before discussing our results, we first briefly summarise 
models which produce 2a galaxies. Bois et al. (2011) used 


idealised numerical simulations to study the mergers of Sb- 
Sc galaxy pairs with dissipation over a range of initial con¬ 
ditions. In the simulations, and in agreement with |Crocker| 
et al. (2009), 2cr galaxies are produced when the merging 


progenitor pair have opposite spin and are coplanar. |Tsatsi| 
et al. (2015) recently showed that major mergers of spi¬ 


ral galaxies on prograde orbits could also produce galaxies 
with 2cr features, although in their simulations, the remnant 
galaxy showed a centrally elevated velocity dispersion-a fea¬ 
ture not normally associated with 2cr galaxies. In all cases, 
galaxies with 2cr features are flattened (e > 0.4) but could be 
either fast or slow central rotators. The extent of the kine¬ 
matically distinct component (KDC) depends on the initial 
orbits and mass ratio of the merging pair. Progenitors on 
anti-parallel orbits produce counter-rotating cores usually 
within 1 Re while those on parallel orbits produce more ex- 
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Figure 10. Line of sight velocity distribution (LOSVD) for GCs 
in radial bins, and by GC subpopulation. Panels from top to bot¬ 
tom show all the GGs, red GGs and blue GGs respectively. Differ¬ 
ent colours correspond to different radial bins. The middle panel 
reveals a kinematic difference between the inner and outer bin 
red GGs. In the bottom panels the blue GGs have “peaked”, sym¬ 
metrical LOSVDs suggesting radially-biased orbits. 


tended decoupled features. The KDC is usually illusory (van 
den Bosch et al.|2008 ), in the sense that it is formed by the 


superposition of a pair of counter-rotating stellar systems. 
The inclination of the merger remnant is important in ob¬ 
serving the 2cr feature, especially for 1:1 mergers. Though 
the idealised binary mergers of Bois et al. (2011) produced 
2cr-like remnants, the simulations were not cosmological. 
They did not include effects from disc regrowth through cold 
accretion flows, bulge growth through minor mergers or late 
major mergers. These could produce remnants that would 
be intermediate between fast and slow rotators. 

Recently, Naab et al. (2014) used cosmological simula¬ 
tions to track the assembly history of present day galaxies 


from 2 : ~2 within the “two-phase” galaxy formation model. 
The only galaxy, out of 44 in the simulations (~ 2%), with 
a clear 2cr velocity dispersion feature (their M0209 in class 
E) is a slow rotator that experienced a late {z < 1) gas-poor 
major merger, and very little in-situ star formation (since 
z ~2). It is therefore dominated by old (~10 Gyr) stars. We 
suspect that the shortfall in the reported number of remnant 
galaxies with a clear 2cr feature could be due to multiple late 
mergers which could conspire to wash away the 2a feature. 

The two simulations described above do not consider 
GCs. There is a dearth of simulation studies for GC kine¬ 
matics in ETGs compared to stellar kinematics. [Bekki et al.| 
(2005) studied the kinematics of GCs in disspationless major 
mergers of MW-like galaxies using numerical simulations; no 
new GCs were formed in these simulations. Even though the 
GCs in the progenitor spirals were not given any initial rota¬ 
tion, they were predicted to have significant rotation at large 
radii in binary major merger remnants. This was interpreted 
as conversion of orbital angular momentum to intrinsic an¬ 
gular momentum. Also, the blue GCs were predicted to have 
a higher central velocity dispersion than the red GCs with 
the rotational dominance parameter (Vm/a) for GCs interior 
to 6 Re greater than that of the galaxy stars within 2 Re. 
Vm is the maximum rotation velocity, so we take the ratio 
VmIcF as an upper limit to our Erot/cr. The velocity disper¬ 
sion is expected to decrease with radius in binary mergers, 
unlike in multiple mergers, where it is expected to have a 
more flattened profile. 

We therefore compare our results with these model 
predictions, in an attempt to unravel the unique assembly 
history of NGC 4473. 


5 DISCUSSION 


It has been shown from stellar kinematics within 1 Re that 
NGC 4473 has two embedded, counter-rotating stellar discs, 
with mass rati o 3:1 ( Cappellari McDermid|2005| Cappel 


lari et al.|20Q7 ). Counter-rotation beyond 60" along the pho¬ 


tometric major axis, as well as multiple kinematic compo¬ 
nents up to 3 Re, have also been reported from the SLUGGS 


survey (Eoster et al. 

2013 

Arnold et al.||2014 Eoster et al. 

2015). Koleva et al. 

(2011 

) showed that NGC 4473 has a 


younger and more metal-rich stellar population in the cen¬ 
tral parts compared to the average for the galaxy. This 


agrees with the results of Kuntschner et al. (2006 2010) and 


Pastorello et al. (2014), who found an extended, metallicity- 


enhanced region along the photometric major axis. These 
studies, along with the central extra-light from |Kormendy| 
et al. (2009) and Dullo & Graham (2013), suggest that a 


gas-rich event can be linked to the complex features we see 
in the kinematics of NGC 4473. 

Our results show that the GC system kinematics of 
NGC 4473 is equally as complex as that of the galaxy stars, 
suggesting co-evolution between the GCs and galaxy stars. 
Within the GC system, the blue and red GCs independently 
show distinct and complex kinematics suggesting different 
formation processes and/or epochs. This again reinforces the 


conclusion that GC colour bimodality is real (Brodie et al. 
2012 Usher et al.||2012| ) reflecting some more fundamental. 
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underlying galaxy property i.e. metallicity, as well as the 


complex formation and assembly history of galaxies (Tonini 


2011) and the radial extent of the counter-rotating region 


2013). The velocity dispersion and the kinematic position 
angle profiles of the red GCs trace those of the galaxy stars 
in the region of overlap, unlike the blue GCs (see Fig. 
an d[^. Similar features are seen in some other galaxies (e.g. 


Schuberth et al. |2010|[Coccato et al.|2013l|Pota et al.|2013 |, 

where the red GCs and galaxy stars are believed to share a 
common formation history. 

Both GC subpopulations however show a large offset 
between their rotation amplitudes and that of the galaxy 
stars. This is odd, as studies of GC system kinematics usu¬ 
ally find a match between the rotation profiles of stars and 
the red GCs (e.g Coccato et al.|2Q13 Pota et al.|2013 ). Here, 
contrary to the norm, we find that the rotation profiles only 
match when all the GCs i.e. the blue and red GC subpopu¬ 
lations are used together in our analysis. Having shown that 
our Hrot profiles are not biased, we suggest that this off¬ 
set is due to the presence of the two counter-rotating, flat¬ 
tened stellar structures at the galaxy centre. When galaxy 
starlight from the radially overlapping stellar components is 
integrated along the line of sight, as in the case of the stel¬ 
lar kinematics from SAURON and SKiMS, the net effect is 
to lower the recovered rotation amplitude, since the compo¬ 
nents are counter-rotating. When we combined all the GCs, 
i.e. without separating them into blue and red GC subpopu¬ 
lations (see Fig.|^, we observe that the rotation amplitude is 
significantly reduced, comparable to what is obtained for the 
galaxy stars. Here, separating the GCs into subpopulations 
by colour also separates them into counter-rotating compo¬ 
nents and reveals the true rotation amplitude. We therefore 
predict that similar offsets would be obtained between the 
rotation profiles of galaxy stars and discrete tracers for all 
2a galaxies, depending on the luminosity contributions from 
the co-spatial stellar components. 

We have also found counter-rotation in the GCs, with 
a sharp transition in both blue and red GCs, at ~ 4 Re- 
This adds to the kinematic complexity of the GC system. 
In the inner ~ 4 Re, blue and red GCs counter-rotate along 
the photometric major axis, while beyond this, the blue 
and red GCs show a sharp ~ 180° kinematic twist, with 
the blue GCs rotating intermediate to the major and mi¬ 
nor axes in the GC system outskirts. This transition ra¬ 
dius and counter-rotation, unambiguously show that the 
galaxy mass assembly occurred in phases from materials 
that had decoupled angular momenta. This is in agreement 
with the “two-phase” galaxy formation paradigm and pre¬ 


dictions from Crocker et al. (2009) and Bois et al. (2011). 


We find further evidence for distinct kinematics on either 
side of this transition radius from the colour profiles of kine¬ 
matic parameters and the LOSVD (see Figs. and 10). We 
note that this strong internal transition is however not seen 
in photometry, underscoring the importance of kinematics 
in obtaining a complete picture of galaxy assembly. 

Gas-rich major mergers have been shown to produce 
sharp radial kinematic transitions, similar to what we have 
observed, between 1 and 5 Re (Hoffman et al. 2010). Ko- 


with coherent GC kinematics we observe here argues against 
a minor merger origin. While counter-rotating stellar com¬ 
ponents can be formed from gas accreted on retrograde or¬ 
bits, this channel is unlikely to be responsible for the large- 
scale counter-rotation observed in the GCs, unless the ac¬ 
cretion event was early (Algorry et al.|2014 Danovich et al 


2014 Lagos et al.||2014 ). Explaining the sharp radial kine¬ 

matic transition and counter-rotation of the blue and red 
GCs in a 2a galaxy in terms of minor mergers is difficult. 
We therefore link the large-scale stellar and GC counter¬ 
rotation to the gas-rich major merger event that formed 
the 2a feature. 

The classic gas-rich major merger scenario for GC for¬ 
mation (Ashman & Zepf 1992) suggests that the merger 


remnant, apart from acquiring GCs from the progenitor spi¬ 


rals, can also form metal-rich stars and GCs (Bournaud 
et al. |2008| [Kmijssen et al. 2011). The group of dynami¬ 
cally “hot”, very red GCs in the inner 4Re radial region 
(see Fig.|^ could have been formed in such a merger event. 
Their rotation direction suggests that they are kinematically 
coupled to the galaxy stars. They could be the GC counter¬ 


part of the younger stars reported in Koleva et al. (2011). 


Measuring their age, metallicity and alpha abundance could 
reveal a younger GC population that matches the ~6 Gyr 


age measured by Koleva et al. (2011) in the galaxy core. 


The central velocity dispersion of the blue GCs is lower 
than that of the red GCs, contrary to the prediction of the 
major merger simulation of Bekki et al. ( 2005| (see middle 
panel in Fig. |^. We do not see signs of increasing rotation 
amplitude or rotation dominance with radius, for both blue 
and red GC subpopulations as predicted. However, our ro¬ 
tation profiles agree with the prediction of fast rotation at 


large radii from Bekki et al. (2005), with comparable Vvot/cF 
for blue and red GCs at all radii. It should, however, be 
noted that the predictions from Bekki et al. ( 2005| are for 
a dry major merger, where no new stars are formed. It has 
been shown from simulations (e.g. Barnes Hernquist|1996 


Jesseit et al.|2007 Hoffman et al.j20T0 ) that the presence of 

gas significantly alters the kinematics of merger remnants, 
especially in the central parts, alongside the formation of 
new GCs. 

In the inner region, where RAkin is aligned with the 
photometric major axis (see Figs. and |^, the kinematic 
misalignment is negligible i.e. ~0°. However, in the outer re¬ 
gion, where the blue GCs dominate, the average RAkin is in¬ 
termediate between the photometric major and minor axes, 
implying a non-zero kinematic misalignment. This suggests 
that the GC system of NGC 4473 is triaxial, closely following 
the galaxy (see also Foster et al.||2013 ). Orbital structures 
from the merging progenitors are predicted to be preserved 
in the galaxy outer regions, where the effect of dissipation 
is negligible. In the galaxy outskirts, the “peaked” shape of 
the LOSVD of the blue GCs suggests that they are pref¬ 
erentially on radial orbits. This agrees with the predicted 
multiple minor merger or accretion origin for the blue GCs 


rmendy et al. (2009) suggested that the counter-rotating ( Cote et al.|1998| Hilker et al.||1999 Dekel et al.|2005 Bekki 


component in NGC 4473 might have formed in a late gas 
accretion or gas-rich minor merger event, due to the small 
region (~ 0.8 Re) associated with the extra-light. However, 


gas-rich minor mergers do not form 2a galaxies (Bois et al. 


et al.||20Q8 ). Incidentally, the same argument can be made 
entirely from the Vrms profiles (see Fig. [^, where the flat 
profile for the blue GCs indicates a multiple minor merger 
origin ( Bekki et al.||2005 Pota et ^|2013[ ). This, therefore. 
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suggests that the outer envelope of the galaxy was assembled 
through minor mergers. 


6 CONCLUSIONS 

Here we present mean velocity and velocity dispersion 2D 
maps from GCs and stars in NGC 4473 to reveal complex 
kinematic features which extend up to 10 We show that 
the double sigma feature of the stellar velocity dispersion 
map extends to the GG system, reaching out to 3 Re, but is 
misaligned with respect to the photometric major axis. 

By fitting an inclined disc model to our GG data, we 
find that the blue and red GGs have different kinematics, 
with a sharp transition in the kinematics at ~ 4 in both 
GG subpopulations. In the inner (< 4 Re) region, the blue 
and red GGs counter-rotate along the photometric major 
axis while in the outer (> 4 Re) region, the blue GGs rotate 
intermediate to the photometric major and minor axes. The 
red GGs in the outer region rotate along the major axis, but 
in opposite sense compared to the inner red GGs. The large 
scale GG counter-rotation, sharp kinematic transition, both 
as a function of galactocentric radius and GG colour, and 
the discovery of a group of centrally located, kinematically 
hot, very red GGs provide additional evidence that NGG 
4473 could have been formed in a gas-rich major merger 
event. In the future, stellar population analysis of this hot 
GG subsample would help determine the merger history of 
the galaxy. We find that the outer region is dominated by 
GGs on radial orbits with a flat Urms profile - which suggests 
that multiple minor mergers may have contributed to the 
build up of the outer regions of NGG 4473. 

We conclude that the GGs in NGG 4473 share a com¬ 
mon assembly history with the stars, based on the similar¬ 
ities in their radially decoupled, complex kinematics. Our 
results differ significantly, in some aspects, from the predic¬ 
tions of GG kinematics in simulated mergers of spiral galax¬ 
ies. We attribute this to the assumptions made in the simu¬ 
lations. A detailed simulation of GG kinematics in galaxies 
formed via gas-rich mergers with dissipation is needed to 
better understand rare but interesting galaxies like NGG 
4473. 
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